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Differential rotation of main-sequence dwarfs and its 
dynamo-efficiency 
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ABSTRACT 

A new version of a numerical model of stellar differential rotation based on mean- 
field hydrodynamics is presented and tested by computing the differential rotation of 
the Sun. The model is then applied to four individual stars including two moderate 
and two fast rotators to reproduce their observed differential rotation quite closely. 
A series of models for rapidly rotating (P rot = 1 day) stars of different masses and 
compositions is generated. The effective temperature is found convenient to parame- 
terize the differential rotation: variations with metallicity, that are quite pronounced 
when the differential rotation is considered as a function of the stellar mass, almost 
disappear in the dependence of differential rotation on temperature. The differential 
rotation increases steadily with surface temperature to exceed the largest differential 
rotation observed to date for the hottest F-stars we considered. This strong differen- 
tial rotation is, however, found not to be efficient for dynamos when the efficiency is 
estimated with the standard Cn -parameter of dynamo models. On the contrary, the 
small differential rotation of M-stars is the most dynamo-efficient. The meridional flow 
near the bottom of the convection zone is not small compared to the flow at the top in 
all our computations. The flow is distributed over the entire convection zone in slow 
rotators but retreats to the convection zone boundaries with increasing rotation rate, 
to consist of two near-boundary jets in rapid rotators. The implications of the change 
of the flow structure for stellar dynamos are briefly discussed. 
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1 INTRODUCTION 

The theory of stellar differential rotation is mainly focused 
on the Sun where helioseismology provides all the possibil- 
ities for a detailed testing of computations. Applications to 
other stars are, however, temp ting in view of the rapid devel- 
opment of asteroseismology (|Christensen-Dalsgaardl |2008| 1 
that can eventually provide data on internal stellar rotation 
(|Suarez et al.ll2010h . A knowledge of differential rotation is 
also considered as a key for stellar dynamos. 

Until recently, precise measurements of differen- 
tial rotation were mainly provided by Doppler imag- 
ing (iDonati fc Collier Cameronl Il997l : IStrassmeieil 12004 



iBarnkeLdJ bd). T l^le^els most suitablTTr 



rapidly rotating stars |Donati 1996ij). Young dwarfs with ro- 
tation periods P rot ~ 1 day present, however, difficulties 
for theory. Rapid rotators have thin boundar y layers at the 
top and bottom of their co nvection zones (jDuimevi fl989: 
iKitchatinov fc Rudigerlll999l ), which are difficult to resolve 
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numerically. Almost no computations of differential rotation 
for rapid rotators were attempted. 

The situation has changed recently. First, mea- 
surements of the differential rotation of two not too 
rapidly rotating, P rot ~ 10 days, main-sequence dwarfs 
were performed using high-precision photometry of the 
asteroseismological MOST mission (|Croll et al ] 120061 : 
IWalker et all 120071 ) . Second, a new mean-field code was 
developed that can resolve boundary layers in rapidly 
rotating stars; this code is used for the first time to 
produce the results of this paper. We first apply it to 
the Sun to check its ability to reproduce helioseismo- 
logical inversions. Then, the differential rotation of the 
MOST-stars (e Eri and a 1 Ceti) are computed and com- 
pared with observations. We also compute the differential 
rotation of AB Dor and LQ Hya which probably have 
been the most frequent observ ational targets for Doppler 
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We further compute the dependence of differential 
rotation on the surface temperature for rapidly rotating 
stars. On a n observation al basi s, the dependence has been 
studied by iBarnes et al.l (|2005l ). The theoretical and ob- 
servational results are quite similar. They both show a 
rapid increase of surface d ifferential ro tation with temper- 
ature. iReiners fc Schmit j (|2002l , l2003h found that the ro- 
tation of F-stars can be strongly non-uniform. Recently 
Ijeffers fc DonatH (|2008l ) observed a large differential rota- 
tion with a pole-equator lap time slightly above 10 days on 
a rapidly rotating GO dwarf. Our computations suggest that 
shallow convection zones of F-stars can possess even stronger 
differential rotation. This raises the question of whether the 
strong rotational shear implies over-normal dynamo activity. 
Our computations suggest a negative answer. The efficiency 
of differential rotation in generating magnetic fields can be 
estimated by the modified magnetic Reynolds number that 
in dynamo theory is conventionally notated as Cn, 



Cn = 



AQH 2 

Vt 



(1) 



IIKrause fc Radlerlll980h , where AS! is the angular velocity 
variation within the convection zone, H is the convection 
zone thickness and ?7 T is the turbulent magnetic diffusiv- 
ity. The Cn parameter ((TJ) is the ratio of the rate Afl of 
magnetic field production by differential rotation to the rate 
rj T H~ 2 of diffusive escape of the field from the convection 
zone. Our computations show that Cn decreases with stellar 
mass. Contrary to intuitive expectations, a small differential 
rotation of M-stars is more efficient in producing magnetic 
fields than the large rotational shear of F-stars. 

Differential rotation models compute the total veloc- 
ity field including the meridional flow. We discuss also the 
meridional flow structure and its variations with the rotation 
rate. 

The rest of the paper is organized as follows. Section 2 
describes our model that is based on the mean-field hydro- 
dynamics (the mathematical formulation is partly shifted to 
the Appendix). Section 3 presents and discusses the results 
and Section 4 summarizes the main findings. 



2 THE MODEL 



The model is based on mean-field hydrodynamics (|Riidigerl 
1989). It computes jointly the differential rotation, merid- 
ional flow, and heat transport in the convection zone of a 
star. The differential rotation of the model results from the 
angular momentum transport by convection (the A-effect) 
and meridional flow. The mod el is close to its former version 
(Kitchatinov fc Rtidigerl [l999) and will be described here 
only briefly. 

In order to compute differential rotation, the model 
needs a knowledge of the structure of a (non-rotating) star 
to specify the basic input parameters such as stellar radius, 
R, luminosity, L, and mass, M. The structure model also 
supplies the density, p c , and temperature, T c , at some small 
depth inside the star, which depth defines the external spher- 
ical boundary (of radius r c ) of the simulation domain. Dis- 
placing the external boundary by 2 — 5 per cent in stellar 
radius below the photosphere helps to avoid problems with 
very sharp near-surface stratification. The reference stratifi- 



cation of the model is adiabatic and spherically symmetric, 
and deviations from the reference atmosphere are computed 
in the model. The deviations are assumed small, so that 
the convection zone stratification should be only slightly su- 
pearadiabatic and the rotation of the star should not be too 
rapid, f2 2 i? 3 (GAf) _1 < 1. Using the values of p c and T c 
as boundary conditions, the equations for adiabatic profiles 
are integrated numerically inwards up to the point where 
the radiative heat flux, 



praa - 



3np 



-VT, 



(2) 



corresponds to the total luminosity, F rad = L(47rr 2 )^ 1 . This 
point is the base of the convection zone. The inner bound- 
ary (n) is placed very slightly (normally by 0.1 per cent 
of the radius) above the base. Opacity n in ([2]) is computed 
using the OPAL opacity tables. Therefore, the input param- 
eters have to include the mass fraction of Hydrogen, X, and 
metallicity, Z. 

So defined, the reference atmosphere helps to specify 
the depth profile of the convective turnover time, 



4c pP i 2 r x 1/3 

3g5F 



(3) 



where c p is the specific heat at constant pressure, g is grav- 
ity, £ — a MLT P(pg)~ 1 is the mixing length, and SF = 
L(4nr 2 )^ 1 — ^? rad is the 'residual' heat flux that convec- 
tion has to transport. The turnover time {J3j enters into the 
key parameter of the Coriolis number, 



fi* = 2rSl. 



(4) 



The differential rotation results fro m the inte r action 
between convection and rotation (|Lebedinskvl Il94ll : 
iTassoul fc Tass oul 2004). The Coriolis number ((4]) measures 
the intensity of the interaction. The value of Q* defines 
whether convective eddies live long enough for the Corio- 
lis force to affect them considerably. 

The angular momentum and heat transport in the con- 
vection zone depend on the Coriolis number. In particular, 
the eddy conductivity tensor, 



that controls the convective heat flux, 
dS_ 

3 



(5) 



(6) 



includes the rotationally induced anisotropy and quench- 
ing via the functions </>(fi*) and </>ii(f2*); explicit expres- 
sions for the quenching functions are given in Kitchatinov, 
Pipin & Riidiger (1994). In (© and (©, S is the specific 
entropy and f2 = fifi -1 is the unit vector along the rota- 
tion axis. Anisotropy of the eddy conductivity (0) means 
that the eddy heat flux is inclined to the radial direction 
even if the entropy varies mainly in radius. As a result, 
the polar regions are warmer than the equator. This 'dif- 
ferential temp erature' is very important for differential ro- 
tation models (|Riidiger et al.| [2005: Mie sch. Brun fc Toomrel 
2006; iBrun fc Rempell 120091 1. It results in deviations of the 
isorotational surfaces from a cylindrical shape. The differen- 
tial temperature on th e Sun has been recently observed by 
iRast. Ortiz fc Meisnerl j2008l V 
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We want to note that the eddy conductivity x T ' s n °t 
prescribed, but expressed in terms of the entropy gradient, 



Xt 



Ti 2 g dS 
' 12c p dr ' 



(7) 



and the same for eddy viscosity. This involves an additional 
nonlincarity in the equations but avoids arbitrary prescrip- 
tions of the diffusivity profiles. The only free parameter of 
our model are c x of ([5|. The former version of the model 
did not use even this parameter, instead assuming that 
c x = 1 (jKitchatinov fc Riidigerll 1999h . In this paper, we put 
c x = 1.5 for closer agreement with helioseismology. 

The model solves the steady equation for the mean ve- 
locity, u, 



(u • V) u + -VP - g = -- V • ipQ) 
P ' P v ' 

together with the entropy equation, 
V ■ (F conv + F rad ) + pTu ■ VS = 0. 



(8) 



(9) 



In ((Hj), Q is the correlation tensor of the fluctuating velocities 



= (?4(T-,£X-(r,t)). 



(10) 



The correlation tensor can be split into a non-diffusive part, 
Q A , representing the A-effect of non-visco us transport o f 
angular momentum by rotating turbulence (tRiidiger]ll989l ). 
and the contribution Q v of eddy viscosities, 



Qij — 



r - — a/" dUk 



(11) 



where AT is the eddy viscosity tensor. The mean flow is as- 
sumed to be axisymmetric about the rotation axis, 



u = ej,r sin 9 SI H — V x 

p \ /' Sll) f 



(12) 



where r, 8, <j> are the usual spherical coordinates, e0 is the 
unit vector in the azimuthal direction, and tp is the stream 
function of the meridional flow. A complete representation 
for the A-effect, eddy viscosities, and two components of the 
motion equation ((8| that provide the equations for angular 
velocity and meridional flow, is given in the Appendix. 

The thermal condition on the top boundary assumes 
black-body radiation of the photosphere. Application of the 
condition on the external boundary is not straightforward, 
however, because of the thin near-surface layer excluded 
from the simulation domain. We assume this layer to be a 
perfect heat exchanger (infinite x T ) so 1 hat the entropy dis- 
turbances at its base and surface are equal. Assuming further 
that the disturbances are produced at constant pressure, we 
find the boundary condition 



F r = 



L 

4nr"£ 



1 + 1 



at r = r e 



(13) 



F r is radial component of the total (convective plus 
radiative) heat flux. The thermal condition at the in- 
ner boundary includes the gravitational darkening effect 
i|R(idiger fc Kukerll2002j ), 



Aivrf 



e + 3 



3cos z 0- 1) at r = n, (14) 



n(n,0)sin 3 d6, 



where Qi is the mean angular velocity and gi is the gravity 
at the inner boundary. 

At both boundaries, stress-free and impenetrable con- 
ditions are applied, 



V> = o, 



r e = at 



n and r 



(15) 



The stress-free condition is a source of certain difficulties. 
The conditions are incompatible with the Taylor-Proudman 
balance of the bulk of the convection zone (the 'thermal 
wind balance', in the geophysical literature). As a result, 
thin layers wher e the balance is violated are formed near 
the boundaries dDurnevlll989l). The boundary layers were 
found in both 3D |Miesch. Brun fc Toomrdl2006h and mean- 
field simulations ( Kitchatinov fc Rudigerlll999l ) of differen- 
tial rotation. The layer thickness is estimated by the Eckman 
depth, D ~ yj f T /(2fi). In the case of rapidly rotating stars, 
the layers can be very thin and difficult to resolve numeri- 
cally (the thickness D decreases with angular velocity faster 
than f2 -1//2 due to the rotational quenching of the eddy vis- 
cosity). 

To resolve the boundary layers, we apply a non-uniform 
grid in radius (zeros of Chebyshev polynomials) with small 
spacing near the boundaries 



1 / / s ( 3 - 3/2 

r j = r <= + n — (r c — n) cos 7r — 

Z \ \ n — Z 



(16) 



2 < j < n - 1, 



ri = n, r n = r c 



n is the total number of grid points. For the latitude depen- 
dencies, Legendre polynomial expansions were applied. This 
leads to the two point boundary value problem in the radius 
for a system of ordinary differential equations. The problem 
was solved numer ically by the standard relaxation method 
l|Press et al.ll 19921 ). 



3 RESULTS AND DISCUSSION 
3.1 Test case: the Sun 

Fig. [1] shows the internal solar rotation computed with 
our model. The figure includes the tachocline region and 
the deeper radiation zone just for completeness of the 
picture. The tachocline was comput ed with a separate 
model (|Riidiger fc Kitchatinov! [2007a) that uses the re- 
sults of the computation of the differential rotation of 
the convection zone as a boundary condition but does 
not influence that computation in any way. The results 
of Fig. [T] are similar to the hel i oseismological rot ation law 
llWilsonTBurtonclav fc TUl 19971 ; ISchou et alii 19981 ) . 

The computed meridional flow is shown in Fig. [5] The 
direction and amplitude of the surface flow are close to ob- 
servations (Komm, Howard & Harvey 1993). Note that the 
flow at the bottom is not small compared to the surface. 
This is a quite general result also found in computations for 
other stars. The stagnation point in Fig.[2]is close to the bot- 
tom so that the downward increase of density does not lead 
to a slow deep circulation. Below the convection zone, the 
meridional flow rapidly decreases with depth. The distance 
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Figure 1. Angular velocity isolines (left) and depth profiles of 
the rotation rate for several latitudes (right) for the model of the 
solar differential rotation. 





0.75 0.80 0.S5 0.90 0.95 
FRACTIONAL RADIUS 



Figure 2. Stream-lines of meridional flow (left) and radial profile 
of meridional velocity for 45° latitude (right) for the solar model. 
Negative velocity means poleward flow. 



of the flow penetration into the t achocline is small com- 
pared to the tachocline thic kness l|Gilman fc Mi esch 2004; 
iKitchatinov fc Rudigerll2006D . 

The meridional flow is closely related to the Taylor- 
Proudman balance. This can be seen from the meridional 
flow equation 



on 2 



g 8S 



VW = sinSr — - JU — 



(17) 



ijKitchatinov fc Riidigeri ri999). In this equation, the left side 
describes the viscous drag due to the meridional flow (its re- 
lation to the eddy viscosity tensor is given by (|A6|) and pip ). 
d/dz = cos6d/dr — r -1 sin Od/dO is the spatial derivative 
along the rotation axis. The two terms in the right side of 
(|17p represent the meridional flow driving by nonconserva- 
tive parts of centrifugal and buoyancy forces, i.e., the cen- 
trifugal and baroclinic drivings of meridional flow, respec- 
tively. The characteristic value of each of these two terms is 
large compared to the left side. Accordingly, the two terms 
nearly balance each other in the bulk of convection zone. 
Consequences of the Taylor-Proudma n balanc e for t he solar 
rotation law have been analyzed by l lDurnevI (|l999T ). Devi- 
ations of isorotational surfaces from a cylindrical shape are 
possible due to the latitudinal inhomogeneity of entropy, 
which in turn can result from the anisotropy of the convec- 
tive heat transport. The Taylor-Proudman balance is cur- 
rently rediscussed in the context of a new hypothesis on the 
coincidence of the isorotational and isentropic surfaces in 
rotating stars (jBalbusll2009l ; iBalbus et al.ll2009l ). 

Fig. [3] shows the contributions of the centrifugal and 
baroclinic terms in (|17[1 , and their sum for the present model 
of the solar rotation. The bulk of the convection zone is 
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Figure 3. Right: Dependencies of the baroclinic (dashed) and 
centrifugal (dashed-dotted) terms of Q17I and their sum (full line) 
on the radius for latitude 45°. The plotted values are normalized 
so that the maximum absolute value equals one. The dotted line 
shows level zero. Left: Isolines of the left side of 1171 . Taylor- 
Proudman balance is fulfilled in the bulk of the convection zone 
but not close to the boundaries. 



very close to Taylor-Proudman balance. There are, however, 
boundary layers where the balance is violated and the layers 
are not very thin (cf. Balbus et al. 2009) . Similar results on 
the Taylor-Proudman balance are provided by 3D numerical 
simulations (Miesch et al. 2006; Brun, Antia, and Chitre 
2010). 

In l|17p. we see that deviations from balance produce a 
meridional fl ow. Observations o f the global meridional flow 
on the Sun l|Komm et all 1 19931 ; IZhao fc Kosovichevl I2004T ) 
indicate that certain deviations from Taylor-Proudman bal- 
ance are present. 



3.2 AfOST-stars 

The main problem with applying the differential rotation 
model to individual stars is to specify the (input) stellar 
parameters. Of the two stars - e Eridani and k Ceti - 
whose differential rotat i on was measured u sing the MOST- 
data jCroll et alj|2006l ; IWalker et al.ll2007l ), e Eri presents 
much less difficult ies because all the required parameters 
were estimated bv lSoderblom fc Dapper] (| 19891 ). 

W e used the EZ code of stellar evolution by IPaxtonl 
(2004) to define the structure of a main-sequence star of 
given mass, age, and composition and infer the input param- 
eters for our simulations from the structure model. Hydrogen 
content was fixed to X = 0.7. The parameters used to model 
the differential rotation of MOST-sta xs are given in Ta bic [[] 
The p arameters of e Eri given by ISoderblom fc Dapper] 
(1989) can be closely reproduced by the structure model 
of O.8M0 star with metallicity Z — 0.01 and an age of 
1 Gyr. The parameters of k 1 Ceti are less certain. Those 
used in differential rotat ion measurements (Rucinsk i et al.l 
12004 IWalker et alj|2007h can be roughly reproduced by the 
structure model of a IMq star with Z = 0.02 at the age of 
about 600 Myr. 

Fig. [4] shows the results of differential rotation simula- 
tion for e Eri. The relative magnitude of the surface differ- 
ential rotation can be estimated with the parameter 



pole 



Our model gives the value of a Dn 
to the observational value of a T ,„ = 



(18) 

= 0.1 27 for e Eri clos e 
0.11 (|Croll et alj|2006h . 



Differential rotation of main-sequence dwarfs 5 



Table 1. Input parameters of the differential rotation models for 
MOST-stars. 



Star 


M/Mq 


R/Rq 


L/L Q 


Z 


Age, Gyr 


Prot 


e Eri 


0.8 


0.724 


0.337 


0.01 


1 


11 


k 1 Ceti 


1.0 


0.907 


0.758 


0.02 


0.6 


9 



P ro t is in days. 





40 60 80 
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Figure 4. Angular velocity isolines (left) and surface profile of the 
rotation rate (right) for the differential rotation model of e Eri- 
dani. 



The agreement for k 1 Ceti is not so close: a DR — 0.130 is 
the computati onal value and = 0.09 is the result of 
measurements ( Wa lker et al.| [2007). The simulated rotation 
laws for both 'moderate rotators' are quite similar. The de- 
pendence of the rotation rate on the latitude in Fig. 2] is 
not as smooth as for the solar model. There is a 'peculiar- 
ity' in the surface profile located around the latitude where 
the angular velocity isoline tangential to the inner boundary 
at the equator arrives at the surface. We always find such 
a peculiarity in rotation laws computed for stars rotating 
considerably faster than the Sun. This means that the often 
used approximation 



fi = n cq (l 



(19) 



may not be very accurate. This peculiarity also means that 
moderate rotators are much closer to the strict Taylor- 
Proudman balance than the Sun. This balance is illustrated 
by Fig. [5] The centrifugal and baroclinic terms in the merid- 
ional flow equation (|17|) strictly balance each other every- 
where except for the thin boundary layers. Violation of this 
balance in the layers excites a meridional flow. Accordingly, 
the meridional flow of Fig. [6] is highly concentrated at the 
boundaries. The bottom flow is not small compared to the 
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Figure 5. Same as in Fig. [3] but for the differential rotation 
model of e Eridani. 
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Figure 6. Simulated meridional flow of e Eridani. Stream lines 
are shown in the left panel and the right panel shows the depth 
profile of the meridional velocity for latitude 45° . 



Table 2. Input parameters of the differential rotation models for 
AB Doradus and LQ Hydrae. 



Star 


M/Mq 


R/Rq 


L/Lq 


Z 


Age 


Prot 


AB Dor 
LQ Hya 


0.9 
0.77 


0.803 
0.698 


0.438 
0.273 


0.02 
0.01 


70 
100 


0.514 
1.6 



Age is given in Myr, P ro t — in days. 

top but the flow in the bulk of the convection zone away from 
the boundary layers is slow. We shall see that the boundary 
layers are even more pronounced in rapid rotators. 

There is an inc reasing belief that meridional flow is 
important for solar llChoudhuri, Schiissler fc Dikpatilll995l ) 
and stellar (jjouve. Brown fc BrurfeOlOl ) dynamos. The flow 
structure prescribed in advection-dominated dynamo mod- 
els is, however, different from the flow predicted by stellar 
circulation models. 



3.3 Rapid rotators 

Simulations of differential rotation were performed for two 
young stars - AB Doradus and LQ Hydrae - that seem 
to be the most frequent observational targets among the 
rapid rotators. Other dwarf stars for which the differen- 
tial rotation was measured by Doppler imaging are ei- 
ther not yet settled on the main sequence or their struc- 
ture parameters are difficult to determine. The parame- 
ters used in differential rotation simulations that also help 
to reproduc e closely the observational structure parameters 
of AB Dor dDonati fc Collier Cameron|ll997| ; lOrtega et~aH 
120071; iGuirado Marti- Vidal fc Marcaide! (20081 ) and LQ Hya 
i Kovari et alj|2004 ) are listed in Table l2l 

Figs. and [5] show the modelled differential rotation 
and meridional flow of AB Dor. The computed differential 
rotation measure a DH = 4.37 x 1 0~ 3 is very close to the ob- 
servational value of 4.5 x 10 -3 (|Donati fc Collier Cameron! 



Il997l ). The peculiarity in the surface profile of rotation rate 
discussed in Section 13.21 is even more pronounced in Fig. [7] 
compared to the moderate rotation case of Fig. 2] The pro- 
file can be only roughly approximated by the cos 2 t9-law of 
()19|) . If the approximation is nevertheless used to describe 
the rotation of the stellar spots, it may lead to a seeming 
variation of the differential rotation with time. The spots 
positioned at different latitudes for different observational 
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g 22.50 




Figure 7. Angular velocity isolines (left) and latitude depen- 
dence of the surface rotation frequency (right) for the differential 
rotation model of AB Dor. 
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Figure 8. Simulated meridional flow of AB Dor. The left panel 
shows the stream lines and the right panel presents the depth 
profile of meridional velocity for the latitude of 45°. Negative 
velocity means poleward flow. 



epochs would lead to different q dr suggesting torsional os- 
cillations even for a steady rotation law. 

Observational estimates of the differential rotation 
of the slower rot ating LQ Hya have a wide spread 
(|Barnes et all 120051 ) . We can only state that the differen- 
tial rotation measure of our model, q dr = 1.28 x 10~ 2 , is 
within the range of observational estimates. 

The meridional flow of Fig.[8]shows an extreme concen- 
tration in the boundary layers. The flow consists of two near- 
boundary jets linked by a very slow circulation in the bulk 
of the convection zone. Such a boundary-layer flow is, prob- 
ably, not important for dynamos. However, the distributed 
(solar-type) flow of Fig. [2] may be significant for magnetic 
field transport. The meridional circulation changes from a 
distributed flow (Fig. [2| to the near-boundary jets (Fig. [5} 
with an increasing rotation rate. This change of meridional 
flow may cause a change in the dynamo regime that may be 
the reason for the two separate branches for fast and slow ro- 
tators in the dependence of t he dynamo-cycle period on t he 
stellar rotation rate found bv lSaar fc Brandenburd l|l999f ). 



3.4 Temperature dependence 

Fig. [9] shows the dependence of the surface differential ro- 
tation on stellar mass computed with our model. The com- 
putations were made for young stars just arrived on the 
main sequence and rotating with a period of 1 day. Models 
were produced for the mass range from 0.4Mq to 1.2M® 
(with O.O5A/0 spacing). These computations cover the sur- 
face temperature range from about 3600 K to 6500 K or 
spectral types from K2 to F6 which roughly corresponds 



o 

TJ 



c 
I 



0.10 7 



0.0 




1 .2 



Figure 9. Surface differential rotation of ZAMS-stais with P ro t = 
1 day as a function of stellar mass. Three lines show the results of 
computations for metallicities Z = 0.01 (dotted), Z = 0.02 (full 
line), and Z = 0.03 (dashed). 



to the range for which iBarnes et al.l l)2005h constructed the 
temperature dependence of the surface differential rotation 
detected by Doppler imaging. The computations were made 
for three metallicity values of Z = 0.01, 0.02, 0.03. For a 
given stellar mass, the results depend on the chemical com- 
position, so that for a mass of 1.2Mq , the surface differential 
rotation differs by a factor of about 10 between the cases of 
Z = 0.01 and Z = 0.03. 

The metallicity dependence almost disappear, however, 
when the differential rotation is plotted as a function of sur- 
face temperature (Fig. I10p . So, the effective temperature is 
indeed convenient for parameterizing th e differential rota- 
tion of young stars |Barnes et al.l I2005T ). An even better 
scaling can be found by using the Coriolis number of @. 
The results for different chemical compositions practically 
coincide on the plot of differential rotation as a function of 
the Coriolis number. As the Coriolis number is not directly 
observable, we shall keep using the surface temperature. 

Figs. [5] and [TD] predict that relatively hot convective 
stars can possess strong differential rotations with pole- 
equator lap times shorter than 10 days. This is larger 
than the stronge s t diff erential rotation observed to date 
ijjeffers fc Donatil I2008T ). The question arises whether a 
strong differential rotation implies over-normal dynamo- 
activity. The results of Fig. QT] suggest a negative answer. 
The figure shows the Cn dynamo-number of ([TJ as a func- 
tion of surface temperature. 

The Cn-parameter was computed for the middle of the 
convection zone using the isotropic part of the eddy mag- 
netic diffusivity r/ ^ = Xt4>(Q*) (that coi ncides with the 
eddy conductivity; iKitchatinov et al.|[l994t ) and latitudinal 
differential rotation (the radial inhomogeneity of rotation 
is relatively small). The estimation assumes that dynamos 
of young star s are distributed over their convection zones 
(|Donatilll999h . 

The Cn of Fig. |TT| declines sharply with tempera- 
ture for F-stars indicating that the strong differential ro- 
tation of these stars is not efficient at producing toroidal 
magnetic fields. This is in agreement with the idea of 
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Figure 10. Same differential rotation as in Fig. [9] but shown 
as a function of surface temperature. The results for different 
metallicities of Z = 0.01 (dotted line), Z = 0.02 (full line), and 
Z = 0.03 (dashed line) now differ by little. 
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Figure 11. Cq dynamo- number of JT} as a function of surface 
temperature. The lines of different styles show the results for dif- 
ferent metallicities, Z = 0.01 (dotted), Z = 0.02 (full line), and 
Z = 0.03 (dashed). 



iDurnev fc Latourl |l978t ) that convective dynamos cease to 
operate at about spectral type F6. The Cn increases steadily 
with decreasing temperature. This is because the convec- 
tion slows down in low mass stars to decrease eddy diffu- 
sion. The decline of magnetic eddy diffusivity overpowers 
the decrease of differential rotation to produce a larger Cn 
in smaller stars. The largest Cn belong to M-dwarfs. This is 
in contrast to the common belief that the small differential 
rotation of M-stars cannot be important for dynamos and 
that the magnetic fields of these stars are generated by the 
a 2 -mechanism. The a 2 dynamos produce nonaxisymmetric 
global fields. However, observations favour an axial symme- 
try o f the global magnetic structure of M-stars (|Donati et al.l 
2006). This may be explained by the effect of differential 
rotation, which is small in low-mass stars, but efficient in 
winding magnetic fields. 

The gravitational darkening of (|14|) is not significant 



for our results. Neglecting the darkening effect reduces the 
computed differential rotation normally by less than 1 per 
cent (by several per cent in extreme cases of rapidly rotating 
F-stars). 



4 SUMMARY 

This paper presents the first results of the new mean-field 
model of stellar diff erential rotation, which impro ves on its 
former formulation (Kitchati nov fc Riid iger 1999), to cover 
the case of rapidly rotating stars with P ro t ~ 1 day. 

The model reproduces very closely helioseismological in- 
versions for the internal solar rotation. The simulated merid- 
ional flow at the bottom of the solar convection zone has an 
amplitude of about 10 ms _1 that is not small compared to 
the surface flow. The near-bottom equator-ward flow can be 
important for the solar dynamo. 

The meridional flow in stars rotating faster than the 
Sun is increasingly concentrated in boundary layers near the 
top and bottom of the convection zone as the rotation rate 
increases. We interpret this boundary-layer structure of the 
meridional flow as an effect of thin boundary layers where 
Taylor-Proudman balance is violated. The change of the 
meridional flow from distributed to boundary-layer structure 
may be the reason for t he change of dynamo regim e between 
slow and fast rotators (|Saar fc Brandenburdfl999l ). 

The differential rotation model was applied to four in- 
dividual stars including two moderate (P ro t ~ 10 days) and 
two fast (Prot ~ 1 day) rotators. In two cases, for which the 
structure parameters of the stars are well known, close agree- 
ment with observations was found. In all cases, the com- 
puted rotation laws were not so close to the sin 2 (latitude)- 
profile of equation ()19|) as it is for the Sun. 

The computations for the rapidly rotating (Prot = 
1 day) ^MS-stars show that the surface temperature, T e g, 
is a convenient parameter for the differential rotation: when 
considered as a function of T e g, the differential rotation loses 
its dependence on the chemical composition of a star that 
otherwise can be quite pronounced. The differential rota- 
tion increases with T e s and the rotation rate difference be- 
tween equator and poles can reach almost 1 radday -1 for 
the hottest F-stars we considered. 

This strong differential rotation is, however, not effi- 
cient for dynamos. The standard Co-parameter of the dy- 
namo models of equation {TJ that measures the efficiency 
of toroidal field production by differential rotation decreases 
with T c ff. Contrary to intuitive expectation, the small dif- 
ferential rotation of M-stars is important for magnetic field 
generation. This may be the reason for the closeness of the 
observed magnetic structure of M-stars to axial symmetry 
(|Donati et alj|2006h . 

As a perspective for future work, theoretical construc- 
tion of the dependence of differential rota tion on stellar age 
and temperature based on girochronology (|Barnes!l2008l ) can 
be pointed out. 



ACKNOWLEDGMENTS 

This work was supported by the Russian Foundation for 
Basic Research (projects 10-02-00148, 10-02-00391). 



8 L. L. Kitchatinov and S. V. Olemskoy 



REFERENCES 

Balbus S.A., 2009, MNRAS, 395, 2056 
Balbus S.A., Bonart J., Latter H.N., Weiss N.O., 2009, MN- 
RAS, 400, 176 
Barnes S.A., 2008, IAUS, 258, 345 

Barnes J.R., Collier Cameron A., Donati J.-F., James D.J., 

Marsden S.C., Petit P., 2005, MNRAS, 357, LI 
Bran A.S., Rempel M., 2009, SSRv, 144, 151 
Bran A.S., Antia H.M., Chitre S.M., 2010, A&A, 510, 33 
Choudhuri A.R., Schiissler M., Dikpati M., 1995, A&A, 
303, L29 

Christensen-Dalsgaard J., 2008, IAUS, 252, 135 
Collier Cameron A., Donati J.-F., 2002, MNRAS, 329, L23 
Croll B., Walker G.A.H., Kuschnig R. et al., 2006, ApJ, 
648, 607 

Donati J.-F., 1996, IAUS, 176, 53 
Donati J.-F., 1999, MNRAS, 302, 457 

Donati J.-F., Collier Cameron A., 1997, MNRAS, 291, 1 

1996, IAUS, 176, 53 
Donati J.-F., Collier Cameron A., Petit P., 2003, MNRAS, 

345, 1187 

Donati J.-F., Forveille T., Collier Cameron A., Barnes J.R., 
Delfosse X., Jardine M.M., Valenti J. A., 2006, Sci, 311, 
633 

Durney B.R., 1989, ApJ, 338, 509 

Durney B.R., 1999, ApJ, 511, 945 

Durney B.R., Latour J., 1978, GApFD, 9, 241 

Garaud P., Ogilvie G.I., Miller N., Stellmach S., 2010, 

arXiv:astro-ph/1004.3239 
Gilman PA., Miesch M.S., 2004, ApJ, 611, 568 
Guirado J.C., Marti- Vidal I., Marcaide J.M., 2008, IAUS, 

248, 496 

Jeffers S.V., Donati J.-F., 2008, MNRAS, 390, 635 
Jouve L., Brown B.P, Bran A.S., 2010, A&A, 509, 32 
Kitchatinov L.L., Riidiger G., 1999, A&A, 344, 911 
Kitchatinov L.L., Riidiger G., 2005, Astron. Nachr., 326, 
379 

Kitchatinov L.L., Riidiger G., 2006, A&A, 453, 329 
Kitchatinov L.L., Pipin V.V., Riidiger G., 1994, Astron. 

Nachr., 315, 157 
Komm R.W., Howard R.F., Harvey J.W., 1993, SoPh, 147, 

207 

Kovari Z., Strassmeier K.G., Granzer T., Weber M., 

Olah K., Rice J.B., 2004, A&A, 417, 1047 
Krause F., Radler K.-H., 1980, Mean-Field Magnetohydro- 

dynamics and Dynamo Theory. Akademieverlag, Berlin 
Lebedinsky A.I., 1941, Astron. Zh., 18, 10 
Ortega V.G., Jilinski E., de La Reza R., Bazzanella B., 

2007, MNRAS, 377, 441 
Miesch M.S., Bran A.S., Toomre J., 2006, ApJ, 641, 618 
Paxton B., 2004, PASP, 116, 699 

Press W.H., Teukolsky S.A., Vetterling W.T., Flannery 
B.P., 1992, Numerical Recipes. Cambridge Univ. Press, 
Cambridge 

Rast M.P., Ortiz A., Meisner R.W., 2008, ApJ, 673, 1209 
Reiners A., Schmitt J.H.M.M., 2002, A&A, 393, L77 
Reiners A., Schmitt J.H.M.M., 2003, A&A, 398, 647 
Rucinski S.M., Walker G.A.H., Matthews J.M. et al., 2004, 

PASP, 116, 1093 
Riidiger G., 1989, Differential Rotation and Stellar Con- 
vection. Gordon & Breach, New York 



Riidiger G., Kitchatinov L.L., 2007, NJPh, 9, 302 
Riidiger G., Kitchatinov L.L., 2007, in Hughes D.W., Ros- 

ner R., Weiss N.O., eds, The Solar Tachocline. Cambridge 

Univ. Press, Cambridge, p. 129 
Riidiger C, Kiiker M., 2002, A&A, 385, 308 
Riidiger G., Egorov P., Kitchatinov L.L., Kiiker M., 2005, 

A&A, 431, 345 
Saar S.H., Brandenburg A., 1999, ApJ, 524, 259 
Schou J., Antia H.M., Basu S. et al., 1998, ApJ, 505, 390 
Soderblom D.R., Dappen W., 1989, ApJ, 342, 945 
Strassmeier K.G., 2004, IAUS, 219, 11 
Suarez J.C., Andrade L., Goupil M.J., Janot-Pacheco E., 

2010, arXiv:astro-ph/1004.0609 
Tassoul J.-L., Tassoul M., 2004, A Concise Hystory of Solar 

and Stellar Physics. Princeton Univ. Press, Princeton, NJ, 

p. 215 

Walker G.A.H., Croll B., Kuschnig R. et al., 2007, ApJ, 
659, 1611 

Wilson PR., Burtonclay D., Li Y., 1997, ApJ, 489, 395 
Zhao J., Kosovichev A.G., 2004, ApJ, 603, 776 



APPENDIX A: MOTION EQUATIONS 
Al Reynolds stress 

The Reynolds stress tensor is related to the fluctuating ve- 
locity correlation of (flO|) and (fTT) . Rij — —pQij. The part 
Q^j of the correlation tensor, which represents the A-effect 
of the non- viscous transport of angular momentum in strat- 
ified rotating fluids, reads 



Qij 



(fiEjkl + TjEikl 

(r ■ n 



(Jo(fi*) + al (n*)) x 
- (Ji(fl*)+a/i(«*)) x 



a 2 



(Al) 



where r is the radial unit vector, a is a parameter of con- 
vection anisotropy (a = 2 in all our computations), and 
7 = Cp/cy. Recent discussions of the A-effect can b e foun d 
in lRiidiger fc Kitchatinov! (1200731 ) and lGaraud et all (|2010h . 
The origi n of the expression (IA1 1 for th e A-effect was dis- 
cussed by iKitchatinov fc Riidigei (|2005l ) where expressions 
for the functions Jo, Ji, lo, h of the Coriolis number £1* 
are also given. 

The viscous part of the Reynolds stress is controlled by 
the viscosity tensor of (|11[) . The viscosity is anisotropic due 
to the rotational influence on turbulent convection, 

Mijki = ^ T (</>i(fi*) (SikSji + 8j k 5a) 

+ 5ik&j£li + Sjk&i&i + SkiCliClj^ 

+ 03(fT) SijSki — <^>4(f2*) SijttkQi 

+ <fc(!r) afijfitfii). (A2) 

The viscosi ty quenching funct i ons, <p n (fl*), n = 1, 5, can 
be found in IKitchatinov et alj (| 19941 ). The eddy viscosity u T 
for a non-rotating fluid is expressed in terms of the entropy 
gradient 
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r£ 2 g dS 
15c p dr 



(A3) 



A2 Angular velocity equation 

The azimuthal component of JSJ gives the continuity equa- 
tion for the angular momentum flux, 

d 



1 9 ( 3 n 

r A p or x 

1 d(r 2 Q)d^ 



sm 2 l 



(sin 9 Qe<t>) 
<9(sin 2 9 fi) dip 



r 2 p 



dr d8 psm 2 9 



i:)r 



+ 

= o, 



(A4) 



where the first and the second lines describe angular mo- 
mentum transport by convection and meridional flow re- 
spectively. On using and (|A1[1 - (]A3|) . the convective 
fluxes of angular momentum can be written as follows 



Qo 



. „ Tt 2 g ds f , 



4>2 (n*) sin 6» ( 



cosy r— - sme- 
ar do J 



ttMLT 

7 



Qr 



sin6>cos6) (JifQ*) + ah(£l*)) 
dQ, 



. „ ri 2 g dS . 



+ cos # I cos # r — — — sin 6^--] 

V ar 96V 

- n ^^MtT^ (j (fi*) + aJo(fi*) 

- cos 2 e (Ji(n*) + o/i(n*)))|. (A5) 



With the angular momentum fluxes (|A5I) . (|A4|) governs the 
angular velocity distribution in the convection zone. 



Q"a, ue and a certain combination of diagonal components 
of Q"j are introduced as new dependent variables when solv- 
ing the meridional flow equation. The expressions for these 
new variables in terms of the stream function can by found 
from IJTTJ, dHJ), and gg) . 



A3 Meridional flow equation 

The equation (|17|l for the meridional flow can be found as 
the azimuthal component of the curled equation l[8]l. The 
left part, 

V{tp) = e^k— - — , (A6) 



drj \p dri 

of this equation describes the viscous drag to the meridional 
flow. In spherical coordinates, (|A6[I reads 



Id f 1 d(pr z Q v re ) 



r dr \ pr 2 dr 
1 d ( 1 <9(sin(9 Q" rf) ) 



r 2 89 \ sin 9 d9 
1 dp 8Q v rr 1 d 2 . v 



rp dr d9 r drd9 
cos 9 d , v 
+ r~s^r9d~r^ Qee ^ Q ^> 

+ ^2§g (Q°o + ~ 2 Qrr) ■ (A7) 

The explicit expression for T>(ip) in terms of the stream func- 
tion is very complicated and never used in practice. Instead, 



